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Abstract. We consider possible routes to superconductivity in hydrated cobaltates 
NaxCo02-yH20 on the basis of the t-J-V model plus phonons on the triangular 
lattice. We studied the stability conditions for the homogeneous Fermi liquid (HFL) 
phase against different broken symmetry phases. Besides the -s/S x -s/S-CDW phase, 
triggered by the nearest-neighbour Coulomb interaction V ^ we have found that the 
HFL is unstable, at very low doping, against a bond-ordered phase due to J. We 
also discuss the occurrence of phase separation at low doping and V . The interplay 
between the electron-phonon interaction and correlations near the ^/3xV3-CDW leads 
to superconductivity in the unconventional next-nearest neighbour /-wave (NNN-/) 
channel with a dome shape for Tc around x ~ 0.35, and with values of a few Kelvin as 
seen in experiments. Near the bond-ordered phase at low doping we found tendencies 
to superconductivity with d-wave symmetry for finite J and x < 0.15. Contact with 
experiments is given along the paper. 
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1. Introduction 

Since superconductivity was discovery in hydrated cobaltates NaxCo02.yH20 for 
X ~ 0.35 and y ~ 1.3 PP an enormous amount of attention was focused on this system, 
in spite of the rather modest critical temperature Tc, that follows a characteristic dome 
shape |2] with maximum value Tc ~ 5K around a doping x ~ 0.35. Indeed, since 
cobaltates may be considered as electron-doped Mott insulators with a layered structure, 
where Co ions are in a low spin state (5* = 1/2) on a triangular lattice, it was expected 
that the resonating valence bond (RVB) scenario proposed long ago by Anderson 3J for 
the cuprates could be clearly realized in the cobaltates. 

The importance of strong electronic correlations, already advanced in view of the 
small bandwidth of the {t2g) levels close to the Fermi energy 01 E], was confirmed by 
recent photoemission studies jHl Ej showing for x = 0.3 a reduction of bandwidth by 
a factor two with respect to the calculated ones |5, . In fact, even taking the whole 
t2g manifold, the bare bandwidth is ~ 1.6 eV while from core level photoemission 
spectroscopy |HI, values of the onsite interaction Udd ~ 3.0 — 5.5 eV were estimated. 
Hence, electronic correlations seem to play an essential role, so that Hubbard or t-J 
models were proposed [0] to describe cobaltates. 

Several unconventional pairing channels and mechanisms were proposed for the 
superconducting state: a) pure t-J model RVB based calculations predict singlet d- 
wave superconductivity CHI ^1 ^] with time-reversal symmetry breaking; b) charge 
fluctuations predict a triplet next-nearest neighbour (NNN) /-wave state jlHl Cll CSj ; 
c) spin-triplet /-wave superconductivity was also proposed based on phenomenology 
and symmetry considerations ^Hl, due to the topology of the Fermi surface ^7], from 
weak-coupling studies of a multiorbital Hubbard model jIHI, and considering spin-orbit 
coupling [1111201121]. 

A number of magnetic resonance [221 12^1 121] and /xSR [221 I2E| experiments are 
consistent with unconventional triplet superconductivity and exclude time-reversal 
symmetry breaking [SHj. Although those results seem at the moment not conclusive, 
with some NMR experiments [2ZI indicating the possibility for singlet s-wave 
superconductivity, specific heat measurements j2H| are consistent with the existence 
of a superconducting gap with nodal lines. Based on the results pointing to triplet 
superconductivity with nodal lines, /-wave symmetry appears as a very promising 
candidate (see also j22j for a more detailed discussion). 

Besides superconductivity, other features of the electronic structure point to the 
proximity of other instabilities. Photoemission spectroscopy indicate the presence of 
a pseudogap of ~ 20 meV with a decrease of the density of states at the Fermi 
energy as the temperature is lowered [201 • Also Raman scattering experiments reported 
recently a pseudogap jHI]. Sidebands in the spectra of Eig phonons suggest that the 
pseudogap arises from a charge ordering instability. Such an instability was only 
observed in superconducting samples. A very recent photoemission experiment |H2| 
was also interpreted in terms of the proximity of the system to a charge order phase. 
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Figure 1. (Colour online) Sketch of the obtained generic phase diagram. The 
figure shows the phases HFL, BOP, VS x V3-CDW and PS. Inside the HFL phase 
superconductivity is possible in the shading regions. Tendencies to superconductivity 
are larger in the regions where the intensity of the shadow is larger. For instance, 
making a cut at a given V along the dashed line we find superconductivity with NNN- 
/ symmetry following a dome shape around x 0.38. For low x superconductivity 
is possible with d-wave symmetry for finite J. With decreasing J the onset of the 
BOP (xc) moves toward a; = and, superconductivity, is at the same time, strongly 
suppressed. The sketch was presented for < a; < 0.5 where our theory predicts, in 
agreement with experiments, a paramagnetic metal. For x > 0.5 the present approach 
does not reproduce the observed Curie- Weiss metal. 



From the possible mechanisms for superconductivity discussed above, b) reconciles the 
possibility of a charge ordering instability with an unconventional superconducting state 
characterized by nodal lines in the order parameter. However, estimates of Tc for a pure 
electronic model lead to extremely low values P^. 

As a natural extension we have recently considered the interplay between phonons 
and electronic correlations jT^j. There, the t-V model on the triangular lattice (where V 
is the Coulomb repulsion between nearest neighbours) was proposed for the electronic 
sub-system. The main effect of V is to bring the system closer to a a/3 x charge- 
density wave (CDW) phase where charge fluctuations are strongly increased. Under 
these circumstances, the electron-phonon (e-ph) interaction vertex is renormalized by 
charge fluctuations leading to superconductivity with NNN-/ unconventional pairing 
symmetry around x ~ 035. 

As mentioned above, originally, several theories were developed on the basis of the 
pure t-J model OEIE]. Moreover, recent NMR experiments [3^ point to enhanced 
low frequency spin fluctuations before superconductivity sets in. Therefore, a more 
general study is needed, where both V as well as the antiferromagnetic spin exchange J 
are considered, in order to clarify the interplay or competition of electronic instabilities 
in different parts of the phase diagram. 

We present here such a study, showing that the region where singlet d-wave pairing 
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due to J dominates is well separated from the one where the NNN f-wave paring appears. 

The presence of J triggers a normal state ^^-independent instability at low doping. 
This so-called bond-order-phase (BOP) is found to be mainly of d-wave character and 
dominated by the exchange interaction J. Phase separation (PS) was also obtained at 
low doping. After identifying the parameter region (J-V-x) where the HFL is stable 
we have studied superconductivity. As in the case with J = 0, unconventional NNN-/ 
pairing near the 73 X V3-CDW ph ase following a dome shape around x ~ 0.35, is 
obtained for expected values of J. The values for Tc obtained in the full model are of 
the order of a few Kelvin like in experiments. In addition, for finite J we have found 
tendencies to d-wave superconductivity in the immediate vicinity of the BOP. Figure Q 
gives a sketch of the phase diagram obtained with our theory. 

The paper is organized as follows. In section II the formalism used is summarized, 
in order to allow for a self-contained presentation. In section III the stability conditions 
for the homogeneous Fermi hquid (HFL) phase are discussed. Section IV is devoted to 
superconductivity and its possible channels. In section V, conclusions and discussions 
are presented. 

2. A large-A^ approach for the t-J-V model and Feynman rules 

Here we describe our treatment of the t-J-V model, given by 

H = -t {ctcja + h.c.) + J Yli^i^j ~ + ^ X] ^i'^j (1) 

where t, J and V are the hopping, the exchange interaction and the Coulomb repulsion, 
respectively, between nearest-neighbour sites denoted by {ij). c|o- (^i<T are the 
fermionic creation and destruction operators of holes, respectively, under the constraint 
that double occupancy is excluded, and rii is the corresponding density operator at site 
i. 

Now, we introduce Hubbard operators [34J which are related with the fermionic 
operators as follows 

Xf = cl n, = (X/Ux/^) 

— ^icJ — '-'i — '-'i 

The five operators X[°'' and Xf^ are boson-like and the four operators and 
X^'^ are fermion-like. The names fermion-like and boson-like come from the fact that 
Hubbard operators do not verify the usual fermionic and bosonic commutations rules. 

In previous papers [33 EH EH EHj we have developed a large- expansion for the 
t-J-V model in the framework of the path integral representation for the Hubbard X 
operators. In this approach the X-operators are treated as fundamental objects without 
any decoupling scheme, and hence, problems that arise in other treatments are avoided, 
like considering fluctuations of gauge fields that appear in the slave boson (SB) approach 
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We start our formalism by extending the Hamiltonian of the t-J-V model, to 
channels for the spin degrees of freedom, and rescaling couplings accordingly. 

+ ^ E ^r^f --"E^r (2) 

{i,j);pp' i,P 

where the spin indices a, a were extended to new indices p and p' running from 1 to A^. 
In order to obtain a finite theory in the A^-infinite limit, we rescaled t, J and V as t/N, 
J/N and V/N, respectively. In (j2)) yU is the chemical potential. 

As shown previously, a path integral can be obtained with an Euclidean Lagrangian 
Le given by 

^E = ^2^ yoo + ^ l3j 

i,p i 

and the following two additional constraints 

^r+E^r-f = o (4) 
p ^ 



and 



ypO yOp' 



xr - = , (5) 

which are required to satisfy the commutation rules of X-operators. While (jH) is the 
completeness condition, originates in the above mentioned requirement to satisfy the 
commutation rules. For a detailed discussion of the constraints, and their relation with 
the commutation rules, we refer to [23110] . In our path integral approach we associate 
Grassmann and usual bosonic variables with Fermi-like and boson-like X-operators, 
respectively. 

We now discuss the main steps needed to introduce a large-X expansion [5^ I37j. 
First, we integrate over the boson variables X^^' using (jSj). The completeness condition 
is enforced by exponentiating and introducing Lagrange multipliers Aj. We write 
the boson fields in terms of static mean-field values, (ro, Aq) and fluctuation fields 
5Aj, as follows, 

= Arro(l + 5Ri) \i = \o + SXi, (6) 

and, we perform the following change of variables for the fermion fields 

Due to (0), and after using (|7j), the exchange interaction contains four fermion 
fields that can be decoupled in terms of the bond variable Aj^ through a Hubbard- 
Stratonovich transformation, where Aij is the field associated with the quantity 
J2p fjpfip/[{^ + 6Ri){l + SRj)]^^"^. We write the Aij fields in term of static mean field 
values and dynamical fluctuations = A{1 + r^ + iA'^), where r^ and correspond to 
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the amplitude and the phase fluctuations of the bond variable, respectively. The index 
rj takes three values associated with the bond directions rji = (1,0), 772 = and 
?73 = (— |, of the triangular lattice. 

Introducing the above change of variables and, after expanding 1/(1 + in powers 
of 6R, we arrive at the following effective Lagrangian: 

Leff = E (f^pfl + fUp) (1 - + ^R') +troY: ifUp + h.c.) 

hp ^ ^ {l,j),p 

-/i E 4/*p(i - + ^R^) + NroY. ^>^^ ^R-^ + E 4/^p(i - ^R^) 



-A' E [i^')" + i^'f] + NrliV - -J) J: 6RM, 



-A Y: ifUp + 4/.p)[l - + SR,) + \sRm, + ^{6R^ + 

-A E [flf,pirH^A'!)[l-\{6R, + 6R,)] + h.c.], (8) 
{hj),p,ri 

where we have changed /i to /i — Aq and dropped constant and linear terms in the fields. 

Looking at the effective Lagrangian (8), the Feynman rules can be obtained as 
usual. The bilinear parts give rise to the propagators and the remaining pieces are 
represented by vertices. 

To leading order in 1/A^, we associate with the A^-component fermion field fp a 
propagator connecting two generic components p and p', 



n 

which is of 0{1) and where Ek is 



Gf){k,u^) = --^ (9) 



Ek = — 2(tro + A) (cos + 2 cos cos —^ky) — /i , (10) 



kx V3, 
— cos — 
2 2 

and k and Vn are the momentum and the fermionic Matsubara frequency of the fermionic 
field, respectively. The fermion variables fip are proportional to the X-operators ((Tj) and 
should not be associated with the spinous from the SB approach. 

The mean field values ro and A must be determined by minimizing the leading 
order theory. From the completeness condition ro is equal to x/2 where x is the 
electron doping away from half-filling. On the other hand, the expression for A is 

A = ^^Ecos(fc>i.(i?fc), (11) 

* k,ri 



where is the Fermi function and Ng is the number of sites in the lattice. For a given 
doping X] jj, and A must be determined self-consistently from {1 — x) = -^J2k'^F{Ek) 
and (HD). 

We associate with the eight component boson field 

6X^ = {6R , 6X, r''^ r''^ r'?^ Ai\ A'^\ A^^), 
the inverse of the propagator, connecting two generic components a and b. 
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a) Propagators and vertices 



pp 

P 

'k,v„ 

p' 



_q,CO„ \j( k,v. 



b) n,,= O + ^ 

d;1= (=====>=====)-' = [n;;T'-n. 



Figure 2. Summary of the Feynman rules, a) Solid lines represent tiie propagator G^"^ 
© . Dashed lines represent the 8x8 boson propagator l(T^ for the 8-component 
field 6X°- . Notice that the component (1,1) of this propagator is directly associated 
with the charge operator. KPJ' (|13() and A^|^ represent the interaction between 
two fermions /p and one or two bosons 8X°' respectively, b) Ilab contributions to 
the irreducible boson self-energy. Double dashed lines correspond to a dressed boson 
propagators. 



^(0)a6(q'^n) 



N 



\ 



where 7^ = {W - 1J){xl2f (cosA;^ 
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(12) 



+ 2cos^cos^/cy) 



and the indices a, h run from 



1 to 8. q and ujn are the momentum and the Bose Matsubara frequency of the boson 
field, respectively. 

The first component bR of the field is connected with charge fluctuations 
i.e. Xf^ = NrQ{l + 6Ri), where X™ is the Hubbard operator associated with the number 
of electrons. 

The non-quadratic terms in (jSI) define three and four leg vertices: 
The three-leg vertex 



-1^ 



^('^n + '^'n) + A* + 2A^cos (A;^ - cos^, 1, 



2A cos(A:,,-^) , -2 A cos (fc,, - ^ ) , 
2 A cos[A;,3-^),2A sin - ^ 
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2 A sinik,,-^] ,2A sin ( - ^ 



(13) 



represents the interaction between two fermions and one boson. 

The four-leg vertex A^l represents the interaction between two fermions and two 
bosons. The only elements different from zero are: 

(In + Q-'r,-, \ % I'v , Qv + l'v 
cos — COS + COS — - — - 




^Trsr = \oi'^n + v'n) + /i + A ^ cos(A;^ - 

1 



^Hr5\ — 2 



and 



A?L. = -A cos(fc,-^^) cos I 5^-' 



Af^'^. = Asin(fc,-^^) cos I 5--'. 



Each vertex conserves momentum and energy and they are of (9(1). In each diagram 
there is a minus sign for each fermion loop and a symmetry factor. Feynman rules are 
given in figure Efe-- 

As usual in a large- iV approach, any physical quantity can be calculated at a given 
order just by counting the powers in 1/A^ of vertices and propagators involved in the 
corresponding diagrams. In the present summary there is no mention of the ghost fields. 
They were treated in previous papers [33 EZI and the only role they play is to cancel 
the infinities given by the two diagrams of figure Eb- 

The exchange interaction J enters fll2|) in two different channels: a) The term 2 J in 
the element (1,1) of is due to the charge-like term, —J/2NJ2{i,j);pp' Xf^Xj^ , of the 



t-J-V model. This term has the same form of the Coulomb term V/NJ2{i,j);pp' Xf^Xj ^ . 
b) The terms 4A^/J in the diagonal of D^J^ are due to the exchange-like term, 

^/2A^E(ij>;pp' XfP'xfP, of the t-J-V model. 

In fll2|) V is only present in the element (1, 1) of D'^q^ and it is multiplied by 
which means that its is strongly screened, at low doping, by correlations. In addition, 
the effect of V is diminished when J is finite. 

The bare boson propagator D{Q)ab (the inverse of (|T2|) ) is 0{1/N). From the Dyson 
equation, D^^ = -0(0)^6 ~'^ab, the dressed components Dab (double dashed line in figure 
I2b) of the boson propagator can be found after the evaluation of the 8x8 boson self- 
energy matrix Ilab. Using the Feynman rules Uab can be evaluated through the diagrams 
of figure |2l3. It results 

n„fe(q, tUn) = -— 2^ ha h : -Tf SalSbl 1" nF{Ek) ,(14) 

where 



h 



a 



1, -2A cos (fc,, + ^) , -2A cos (a:,, + ^ 



+ ^) , 2A sin [k,, + ^) , 2A sin [k,,, + ^) , 2A sin [k,, + ^ 
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and Ek = —2t (x/2) Y^ri cos(A;^). 

The component (1, 1) of the dressed boson propagator Dab (also caUed Drr) is 
related with the charge-charge correlation function Xij- K can be written as [SHIISI 

4M = ^E(^^^rW(o)). (15) 

pq 

and the completeness condition and the relation between and 6Ri, in Fourier space, 
X%q,u;) = -N(^^yDnR{q,uj). (16) 

In EH E] it was pointed out that in C(l) the charge-charge correlation function 
shows the presence of collective peaks above the particle-hole continuum. 

Finally, one remark is in order at this point. From the completeness condition 
we can see that the charge operator is of 0{N), while the operators X'^^ are of 
(9(1). This fact will have the physical consequence that the 1/A^ approach weakens the 
effective spin interactions compared to the one related to the charge degrees of freedom. 
Another consequence of this result is the absence, in 0{1), of collective excitations 
(like magnons) in the spin susceptibility. The spin-spin correlation function is then a 
Pauli like electronic bubble with renormalized band due to correlations |^ . While 
there are collective effects in the charge sector in 0{1), those appear in 0{1/N) in the 
spin sector. However, while focusing on superconductivity, that in cobaltates occurs at 
relatively large doping where the system behaves as a paramagnetic metal, this fact is 
not relevant. 



3. Phase diagram: Instabilities of the homogeneous Fermi liquid 

Before considering the possible electronic instabilities, we recall that in leading order 
we have free fermions with an electronic band (jlOj) . renormalized by Coulomb 
interactions. From this electronic dispersion we obtain a large Fermi surface (FS) 
enclosing the F point. First principles calculations [H] predict, apart from this FS, 
the existence of small pockets near K points. However, it is important to notice 
that recent ARPES experiments [Hj do not show the presence of pockets. Invoking 
electronic correlations, a theoretical explanation for the absence of pockets was given 
using LDA+U jl2] and a strong coupling mean field approach [IH]. These results give 
an additional support for considering cobaltates as strongly correlated systems. 

Instabilities of the C(l) HFL phase are studied by analyzing the zeros of Det D^i^ = 
Det {D^Q^^fj — Hab) ■ Expanding the determinant by minors along the first row, it can be 
written in the static limit {un = 0) as follows, 

Det = fi^, J,q)V + g{x, J, q) (17) 

where / and g are long algebraic expressions and they were computed numerically in 
order to study the instabilities. The fact that the Hq^'s are V independent functions 
was used in fITTj) . 

The system presents two kind of instabilities: 
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Figure 3. F-dependent phase diagram for J = 0.2. Vci marks de border between 
the HFL and the \/3 x ^/3-CDW. 142 marks de border between the HFL and PS. The 
vertical dashed line separates the HFL from the BOP (see subsection 3.2^. 



a) V^-dependent instabilities: They occur when / 7^ 0. From (fT7|) these instabihties 
take place for given x, J and q, when the Coulomb potential is V = Vc = —g/f- 

b) ^-independent instabilities: They occur when the two functions / and g are zero 
simultaneously. 

Considering the bare hopping t ~ 150 meV jl] and Udd ~ 3.0 — 5.5 eV ^ we obtain, 
in agreement with other estimates [HE] {J = At'^/U) J ~ 0.1 — 0.2 in units of t, which 
in the following is considered to be 1. In what follows we choose J = 0.2 which can be 
seen as an upper bound value. 



3.1. Charge density wave instabilities 

In this section the stability of the HFL phase will be studied as a function of V, since it 
corresponds to type a) above, i.e. to the V^-dependent kind. Two instabilities are found 
with critical values Vd and Vc2. In figure El we show the phase diagram in the Vc — x 
plane for J = 0.2. Regions corresponding to HFL, a/3 x v^-CDW {V > Vd) and phase 
separation (PS) (V < Vc2) are identified. 

As usual, PS takes place at q = 0. For a given J, the system shows PS for low 
X and V below the corresponding line marked with Vc2. The PS region increases with 
increasing J. 

When V > Vd the system enters in a CDW state. The divergence at the onset of 
the CDW is at q = Q = (4/37r,0) corresponding to a ^3 x VS-CDW [HIISl- The 
Vc — x line separating the HFL from the a/3 x ^73 CDW phase has a parabola-like shape 
with a minimum closer to the doping x ~ 0.35 where superconductivity takes place, 
with maximum Tc, in cobaltates. 

The critical Coulomb repulsion Vd increases with increasing J. As can be seen 
from the element (1, 1) of DTq}^^ in |T2|l . the effect of V is diminished by the presence 
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Figure 4. ^/-independent instability. For a given J the HFL is unstable for x < Xc 
against the BOP. Inset: vs J. 

of J. For instance, for x = 1/3, Vd is 1 and 1.13 for J = 0. and J = 0.2, respectively. 
For J = our phase diagram agrees with the obtained one in (see figure 2 of that 
paper). The eigenvector corresponding to the zero eigenvalue of D~i^ is mainly of the 
form (1, 0, 0, 0, 0, 0, 0, 0) and therefore the instability is concentrated in the charge sector. 

The vertical dashed line separating the HFL from the BOP will be discussed in the 
next subsection. 

3.2. Bond- order phase 

This instability corresponds to the V-independent kind. In figure ID we show the phase 
diagram of the model considering only this instability. For a given J, below a critical 
doping Xc, the HFL is unstable. As can be seen, the instability is strongly J dependent. 
When J 0, Xc showing that J is the main source for the instability. For instance, 
for J = 0.2, Xc ~ 0.066 which corresponds to the vertical dashed line in figure El In 
the inset of figure 0] we show the J dependence of the critical momentum qc where the 
instability takes place. It is of the form qc = {qx,0) and leads to an inconmensurate 
instability. 

The dominant symmetry involved in the instability is given by the eigenvector 
corresponding to the zero eigenvalue, that in this case has the form ~ 
(0, 0, 0, -V2/2, y/2/2, 0, 0, 0), and hence, it is in a different sector than the V^-dependent 
instabilities. It means that, for a given J, the amplitude variables r^^ and r^^ of A'' 
(the fourth and fifth components of SX'^) are frozen at qc for x < Xc- The term 
~^J2{i,j),p,r){fipfjp + fjpfip)^i^ ill the fourth line of our effective Lagrangian (jH)) (from 
which the components 3, 4 and 5 of the vertex A^^' (fT3|) are obtained) is, in g-space, of 
the form 2A J2kqp,ri cos(/i;^ — g^/2)r^/t^_,_^ p/fc^p. When the variables r''^ and r^'-^ are frozen, 
a new hopping-like term of the form [cos{qcx / 4:)dxy{k) + sm{qcx/4:)px{k)] is generated in 
the Hamiltonian. px and dxy are harmonics of the triangular lattice jHj. Hence, the 
instability can be purely of c? or p character depending on whether it takes place at 



Electron-phonon interaction dressed by electronic correlations 



12 



Qc = (0, 0) or qc = {2tt, 0) respectively. In-between, both symmetries are mixed. For 
instance, for J = 0.2 we have qc ~ (0.8, 0)7r which means that the instabihty is ~ 65% d 
and ~ 35% p. Therefore, it is mainly of rf-wave character at the onset of the instability. 
Since in the triangular lattice rf-wave symmetry is two-fold degenerate, the new phase 
can be a combination of both. 

In the case of the t-J model on the square lattice, these kind of purely electronic 
instabilities were studied using SB jH], Bayn-Kadanoff funtional theory |^ and path 
integral large-N approach [SZj. In this case two regimes were obtained: a) For J < 0.5, 
at low doping the system shows an instability whose eigenvector is mainly confined to 
the sector corresponding to phase fluctuations of the bond variables. Therefore, the 
new phase has a complex order parameter and corresponds to the well known flux phase 
(FP). b) For J > 0.5, at low doping the eigenvector of the instability is confined to the 
sector corresponding to the amplitude variables r^, and hence, corresponds to the BOP. 
This is in fact the order we found in the triangular lattice. From the discussion above, 
and the fact that phase separation is also found on the triangular lattice, it seems that 
instabilities expected on the square lattice at high values of J already appear at low 
values on the triangular lattice, a feature that may be due to the larger coordination 
number. We would like to point out that to our knowledge, this kind of analysis has 
not been done before for the t-J model on the triangular lattice. 

Finally, a word of caution about the BOP is in place here. Our figure ID shows 
that at zero doping, the BOP would set in for J finite. However, numerical results 
|ini 1^ indicate that the antiferromagnetic Heisenberg model on a triangular lattice 
displays long-range Neel order for spin S = 1/2. The same discrepancy appears on the 
square lattice with the techniques metioned above. While due to its bipartite nature, the 
antiferromagnetic order is expected to be rather robust on a square lattice, geometric 
frustration on the triangular lattice should render this state much more fragile, such 
that upon doping, RVB-like scenarios appear even more probable in the present 
case. Therefore, since our large-A^ approach shows phases like the FP and BOP at 
low doping on a square lattice jHZ|, as other mean field aproaches (see above), that are 
considered as serious candidates for underdoped cuprates, we expect the results on the 
triagular lattice to be even more trustworthy. 

4. Superconductivity 

Having studied the stability conditions for the HFL under the influence of V and J, we 
consider in this section possible superconducting states. 

In ^H] we have proposed that the interplay between electronic correlations and e-ph 
interaction is relevant for describing superconductivity in cobaltates. To this aim we 
consider the additional electron-phonon Hamiltonian Hp^ + H^-ph where 




(18) 
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Figure 5. a) Total effective paring Vj^j as the sum of a pure electronic mediated 
V*yj^~^ and a phonon mediated y^Jf^ interactions. Solid, double dashed, and dotted 
lines are the propagators for fermions, I?a6, and phonons respectively. In V^J^^ ^ the 
bare e-ph vertex, g (solid circle), is renormalized by electronic correlations as showed 
in b). The last diagram contains a four leg vertex proportional to g which is generated 
when our X-operator approach is applied (see l|221>)- 



and 

^e-p. = ^ E («1 + «0 • (19) 



In order to obtain a finite theory in the ^ oo we have rescaled the e-ph couphng g 
to gl^. 

From Hph ffTS|l we have a free phonon propagator 

of = ^ , (20) 
which is of OiV)- Using the constraint (jS)), the expression (fTUj) for He-ph reads 

where the relations between X'^'^ and fp and between and 5R were also used. Up 
to 0{1/N) only the first term in the expansion of 1/(1 + 5R) is necessary, leading to 

H^.pu = -j= E(4 + a.)(l - 5R,)flf,p . (22) 

V -1 V ip 

Next, we discuss the effective interaction between carriers. Fluctuations in 0{1/N) 
give rise to interactions among fermions. In 0{1/N) there are two contributions to the 
total pairing effective potential T^/j shown in figure The first diagram of figure 

'eff 

in the fermionic propagator (jU)) are not necessary for calculating V^jf in 0{\/N). 

Due to the rescaling of the e-ph interaction, g, superconductivity from phonons 
also appears in 0{1/N) and therefore, it can be treated on an equal footing to 
superconductivity in the pure t-J-V model. 



shows the pairing effective potential, K// ^ , for the pure t-J-V model while the second 
diagram shows the e-ph pairing potential V^ff^. Notice that 1/A^ self-energy corrections 
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The new contribution to the e-ph pairing potential is the vertex (dark triangle) in 
the second diagram of figure which represents the e-ph interaction renormalized by 
electronic correlations. The diagram of figure Eb shows that the renormalization of the 
bare vertex g is due to the electronic correlations of the pure t-J-V model, which will 
be the main contribution to our results. Notice that the e-ph vertex is not renormalized 
by the e-ph interaction. Hence, we assume, as in usual metals, that Migdal theorem is 
valid. 

Using our Feynman rules the renormalized e-ph vertex 7 is 



7(g, k,k ,u}n, i^n, '^n)=9N \ - + 2A— cos cos A; + - — 

X D^\q,co^)A^^ , (23) 

where A^''^ is given by (fT^ . 

In J = case, the renormalized e-ph vertex can be written as 

7(g, k', k",uj^, v'^, vl) = Ng^ | '-{u'^ + K) + Dnji{q,uJn) + I)A/?(g,^n)} , (24) 

where the relevant contribution comes from the charge-charge correlation Dfui fjl6|) . 

Vertex corrections obtained by us are similar to the early calculation of jlHl which were 

used before for studying transport |49^ and isotope effect in cuprates [HU] . 

In weak coupling we can evaluate the effective interactions on the FS, i.e. for 
= ^'n = ~ ^ ^^"^ momcutum q = k — k' with k y k' on the FS. Then, 

the total pairing effective potential is 

\//^/k, k') = K/-/"''(k, k') + K7/'(k, k') , (25) 
where, using the Feynman rules, 

V;7/-^(k, k') = Af D'^^(k - kOAf , (26) 

and 

K7/^k,k') = -^[7*(k,k')f , (27) 

where A = ^^(0) is the bare dimensionless e-ph coupling and A^(0) is the bare 
electronic density of states. In (jTfjl ^* = g. 

In the following we choose A = 0.4. To our knowledge the value of A is not known 
for cobaltates yet. However, recent experiments suggest that it is nonneghgeable jSI]. 
A = 0.4 is of the order of recent estimates 



Without correlations, 7* = 1 and VeJl^ is the usual pairing potential used in BCS 
theory which in conventional metals leads to superconductivity in the isotropic s-wave 
channel. However, in such a case, the characteristic dome shape observed in cobaltates 
would not be expected, because there is no reason for a strongly doping dependent bare 
e-ph coupling A. Recently, Yada and Kontani [22], using a. d — p model for NaCo02, 
found evidence for phonon mediated superconductivity in the s-wave channel, so that. 
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as they pointed out, other effects, such as electronic correlations, are necessary in order 
to stabilize an anisotropic pairing. 

We use the effective potentials to compute the dimensionless effective couplings in 
the different pairing channels or irreducible representations of the order parameter. The 
critical temperatures Tc, can be then estimated from: Td = 1.13ct;cexp(l/Aj), where Uc 
is a suitable cutoff frequency. 

The effective couplings Aj with different symmetries are defined as jHUj : 

. ^ 1 J{dk/\v^\)J{dk'/M)g,,{k')V,ff{k',k)g,{k) 
' (27r)2 J{dk/\v^\)g,{k)^ ^ ^ 

where the functions gi{k), encode the different pairing symmetries (see Table I of ^1] 
for the triangular lattice), and fk are the quasiparticle velocities at the Fermi surface. 
The integrations are restricted to the Fermi surface. Aj measures the strength of the 
interaction between electrons at the Fermi surface in a given symmetry channel i. If 
Aj > 0, electrons are repelled. Hence, superconductivity is only possible when Aj < 0. 

Figure ini shows results for the most relevant symmetry channels as s, d and NNN-/ 
for J = 0.2 and V" ~ 1 close to the ^3 x VS-CDW (see figure ^ . The curves were 
cut at X ~ 0.066 where the BOP instability takes place. Figure IHt corresponds to the 
pure t-J-V model. As expected A*""'"^ > (dashed hne) hence, electrons are repelled 
indicating no tendencies to superconductivity with s-symmetry. 

In the NNN-/ channel (solid line) we have obtained small negative values for 
AiviViv-/ with a shallow minimum around x ~ 0.38 suggesting the possibility for 
superconductivity. However, in the most favorable case X'^^j^^^j ~ —0.04 implying 
that an unrealistic cut-off Uc is necessary in order to obtain a value of a few Kelvin for 
Tc. This feature remains valid even closer to the x ^/S-CDW phase. 

The dotted-dashed line shows results for A^""^"^. Negative superconducting 
couplings are found at small doping, x < 0.15, where A^""^"^ can take robust values 
~ —0.3, indicating that superconductivity in the d-wave channel may be expected at 
such low doping. Our calculations show that A^^"^"^ is independent of V. On the 
other hand, A^~'^~^ is strongly J dependent, vanishing fast when J — > 0. Comparing 
figures El and IHl it can be seen that d-wave superconductivity occurs near the onset of 
BOP which is also of d-wave character and triggered by J. For the doping range where 
superconductivity takes place in cobaltates, x ~ 0.35, there is no indication for rf-wave 
pairing. 

Figure Eb shows results for the e-ph case. A^"^^ (dashed line) would suggest that 
superconductivity could be expected in the s-wave channel around x ~ 0.38 following 
a dome shape as in cobaltates. However, in a Gutzwiller description the s-wave order 
parameter would be exactly zero. In contrast to that the enforcement of the large-iV 
non-double occupancy constraint (^, namely, that only N/2 out of the total states at 
a given site can be occupied at the same time, makes s-wave superconductivity possible 
in spite of strong correlations (see ^| for discussions). 
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^N'NN-f (solid line) shows robust negative superconducting couplings following a 
dome shape around x ~ 0.38. Notice that X'^-NN-f ~ —0.16 which is a factor 4 larger 
than for the pure t-J-V model. NNN-/ superconductivity is not very sensitive to the 
value of J and the only requirement is that the system must be close to a charge 
instability (in this case a x a/3-CDW). 

The dotted-dashed line in figure IHb shows the projection of the e-ph coupling on 
the ci-wave channel, A^"*^^. In the uncorrelated case A^"*^^ is exactly zero, but electronic 
correlations are responsible for the weak modulation with doping of A^~^^. For low 
doping, A^"*^^ can take small but negative values of the order of A^~^'^ ~ —0.02 (not 
well appreciated in the scale of figure EI). These small and negative values for A^"^'' 
increse with J. An appreciable strength of A^~^'* can only be obtained for unrealistic 
values of J. As we mentioned above the renormalized e-ph vertex is mainly dominated 
by charge fluctuations and then, NNN-/ pairing proves correlation effects better than 
d-wave pairing. 

Figure IHt shows results for the total coupling Af = A*~"^~^ + Xl~^^. The solid 
line i^JiNN-f) shows the strongest tendencies to superconductivity. The dome shape 
around x ~ 0.38 is clearly present. While X^^^lv^f is small but negative, A^^^_j 
increases slightly with respect to A^^^_j. 

The situation is different for the s-wave channel; as A*""^"^ > 0, Aj will be smaller 
than Xl~P^. For instance, A^ ~ —0.05 at a; ~ 0.38, as a result of strong correlations. 

Finally, as the contribution A^~^'^ (dotted dashed line in panel b) is very small, the 
total effective coupling A J is close to A|^~"^~^. In the triangular lattice the ci-wave channel 
is degenerate. Our calculation determines the leading symmetry of the superconducting 
order parameter but not its value. However, our results are consistent with previous 
mean field studies [TUl [T2] where superconductivity was found with di + id2 symmetry. 
We would like to remark that these mean-field studies assume the pure J term as 
the effective interaction. In our calculation we have included fluctuations through the 
infinite series of bubbles in the evaluation of the propagator Dah- 

We conclude that the t-J-V model alone would support superconductivity in the 
d-wave channel at low doping [x < 0.15 for J = 0.2). To our knowledge there are 
at present no data for such doping levels concerning superconductivity. On the other 
hand, while e-ph interaction essentially introduces no changes in the d-wave channel, 
it effectively couples to charge fluctuations close to a charge ordering instability. Near 
charge-order, charge-fluctuations renormalize the e-ph effective interaction in such a way 
that superconductivity with triplet NNN-/ symmetry is favored around x ~ 0.38. 

In order to see the influence of the proximity to the charge order on 
superconductivity, in figure d we show, for x = 0.38, the values of X^^^lv^f (dotted- 
dashed line), (dashed line) and A^^^_j (solid line) as a function of V . For the 
pure t-J-V model the values of A^^^^j are very small even very close to the \/3 x \/3- 
CDW. In contrast, when e-ph interaction is included, there is a large window for the 
parameters for which superconductivity may be possible. As A^jvAf-/ increases with 
V , Tc will also increase. It is known from experiments that increases with water 
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Figure 6. Superconducting couplings in the s (dashed line), d (dotted dashed 
line) and NNN-/ (solid line) channels for a) the pure t-J-V model, b) the e-ph model 
and c) the total case. Results are for J — 0.2, and V = 0.914 which locates the system 
near the V3 x \/3-CDW phase. The bare e-ph superconducting coupling was chosen 
to be A = 0.4. 

inclusion supporting the view that the increasing of V mimics the increasing of 
water content. For small V , A^^^_j is very small in agreement with the non-existence 
of superconductivity in nonhydrated samples. 

Although at this point a quantitative comparison with experiments is beyond the 
scope of our analysis, it is still of interest to make a qualitative estimate of T^. Since 
NNN-/ superconductivity is in our case mainly mediated by phonons, we consider as 
a relevant cut-off corresponding energy scale. Recent first principles lattice dynamics 
calculations show the existence of optical phonons as high as 75 meV and for simple 
estimates we consider oje = 40 meV jSH]- Using the values for A^^^_y (solid line in 
panel c of figure IHl) we show, in the inset of figure d results for T^. Our crude estimate 
gives Tc ~ 2K which is of the order of the experimental value ~ bK. A priori, this 
result may appear trivial because the e-ph interaction is certainly efficient for leading to 
superconductivity with in the scale of a few Kelvin in usual metals. However, there 
are two features in our results that would be absent when correlations are left aside. 
Firt, the pairing channel is an unconventional one, a result that would not be possible 
considering e-ph interaction alone. Second, a dome shape is obtained for Tc, again a 
fact that would be missing by considering e-ph interaction alone. 

5. Conclusions and discussions 

We have proposed that the t-J-V model plus phonons, on the triangular lattice, may 
describe superconductivity in hydrated cobaltates. 

Before studying superconductivity we have presented the phase diagram of the 
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Figure 7. Superconducting coupling in the NNN-/-channcl as a function of V 
approaching the critical value Vc for J = 0.2. The figure is presented for x = 0.38 where 
the largest superconducting coupling was obtained. Results for the pure t-J-V model 
(dotted-dashed line), pure e-ph model (dashed line) and the total case (solid line) are 
presented showing that, when phonons arc included, a large window near the V3 X V3- 
CDW phase exists where the system presents indications for superconductivity. Inset: 
Tc vs doping for the NNN-/ channel. The BCS formula was used for the estimation of 
Tc (see text). 

model. Two types of instabilities were found: a) ^-dependent instabilities. These 
instabilities are dominated by the short range Coulomb repulsion V. For V larger than 
a critical value Vc the system enters a \/3 x a/3-CDW phase. PS was also obtained for 
small X and V. b) K-independcnt instability. This instability is dominated by J and is 
common to the pure t-J model. For a given J, the system is unstable for doping smaller 
than a critical one Xc- This new phase occurs at an incommensurate momentum Qc and 
it is called BOP. It was found that — > when J — > 0. 

These phases delimit the region of parameters (x, J and V) where the HFL is 
stable. Paring was calculated in this last region. 

Near the VS x -\/3-CDW state the e-ph vertex is renormahzed by electronic 
correlations developing an anisotropy in /c-space due to the coupling with charge- 
fluctuations. This anisotropy favors superconductivity with NNN-/ symmetry when 
the renormahzed vertex is used for calculating phonon mediated pairing. Besides the 
possibility of anomalous pairing, the model shows possible superconductivity following 
a dome shape around x ~ 0.35 with values for T^. of the order of a few Kelvin. We have 
found that the above results are robust against the value of J around x ~ 0.35. 

In addition to NNN-/ superconductivity for doping x ~ 0.35, for finite J we have 
found the possibility of d-wave pairing at low doping. For instance, for J = 0.2, on 
the high side for possible values of J, d-wave pairing can be expected for x < 0.15. 
In contrast to NNN-/ pairing, d-wave pairing is not affected by V suggesting that 
superconductivity in that channel could exist even without hidratation. To our 
knowledge superconductivity was not investigated for samples at such small doping 
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levels. 

It was recently proposed [52], and studied experimentally [53], that hydration causes 
the electronic structure to be more two dimensional. Notice that in ^2] it was shown 
that Tc decreases with decreasing lattice parameter c. We think that due to this effect 
the Coulomb repulsion V may be less screened when the system is hydrated. If V 
is small when the system is not hydrated, phonons will favor superconductivity only 
in the s-wave channel. However, the strong repulsion in this channel from the t-J-V 
model (figure lEK) suppresses pairing from phonons. This may be the reason for the 
nonexistence of superconductivity in non-hydrated cobaltates. 

As NNN-/ superconductivity has a large contribution from phonons we expect a 
large isotope coefficient. Our theory also predicts a rather constant isotope coefficient 
along the dome in contrast to the strong doping dependent isotope coefficient in cuprates 
[57j . To our knowledge, isotope effect experiments are still not available for cobaltates. 
We expect that the improvement in single crystals preparation |5H] will be useful for 
isotope experiments in the near future. We consider this experiment as a strong testing 
for our approach. 

We would like to remark about some analogies between cobaltates and organic 
materials. As in organic materials [SH], optical conductivity experiments in cobaltates 
[UU] show the presence of low energy features which can be associated with the proximity 
of the system to the charge-order. Merino and McKenzie [HI] pointed out that the 
proximity to the charge-order is relevant for superconductivity with anomalous paring 
in organic systems. As our scenario predicts a \/3 x \/3-CDW state for V > Vc it will 
be interesting to see if a further inclusion of water can trigger the a/3 x a/3-CDW phase 
or, at least, if low energy optical features are reinforced with hydration. 

Recent reports discuss the possibility for singlet s-wave superconductivity and 
the coexistence of s-wave and unconventional pairing 62 . Such a situation could 
be reached in our case by increasing the bare A from A = 0.4, since then, the total 
superconducting couplings Aj and A^^^_j become more attractive and, for A > 1 
both symmetries are nearly degenerated. However, lacking detailed information about 
the e-ph coupling, we take a cautious value for A, that is already sufficient to trigger 
superconductivity, and as shown above, of unconventional type. 
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